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Application of a Linearized Unsteady Aerodynamic 
Analysis to Standard Cascade Configurations 


SUMMARY 


A linearized potent ial- flow analysis, which accounts for the effects of 
nonuniform steady flow phenomena on the linearized unsteady aerodynamic response 
to prescribed blade motions, has been applied to five cascade configurations. 
These include the first, fifth, eighth and ninth standard configurations proposed 
as a result of the Second International Symposium on Aeroelasticity in Turbo- 
machines and a NASA Lewis flutter cascade. Selected results from this study, 
including comparisons between analytical predictions and the experimental 
measurements submitted for three of the foregoing configurations, are described 
in this report. The correlation between theory and experiment for the first 
standard configuration (a compressor cascade operating at low Mach number and 
frequency) is quite good. Moreover, the predictions and measurements for the 
NASA Lewis cascade of symmetric biconvex airfoils show good qualitative agree- 
ment. However, wide discrepancies exist between the theoretical predictions and 
the experimental measurements for the fifth standard configuration (a subsonic 
transonic fan tip cascade). These can be partially attributed to conditions 
being imposed in the experiment which differ from those commonly used in unsteady 
aerodynamic analyses. 
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INTRODUCTION 


In the present study the linearized pot ent ial- flow analysis, described 
briefly in the next section and in detail in Refs. 1 through 4, has been used to 
determine unsteady response predictions for five cascade configurations. In 
particular, detailed unsteady response predictions have been determined for the 
first, fifth, eighth, and ninth standard configurations suggested by T. Fransson 
and P. Suter (Ref. 5) for theoretical and experimental investigations on turbo- 
machine cascades. Similar theoretical results have also been determined for a 
NASA Lewis flutter cascade in support of the unsteady cascade experiments 
conducted at NASA Lewis Research Center (Ref. 6). In this report we summarize 
briefly our contribution to the establishment of the theoretical and experimental 
data base for the standard cascade configurations of Ref. 5. In addition, a 
detailed presentation of our theoretical results for the NASA Lewis cascade will 
be presented along with comparisons between these theoretical predictions and the 
measurements reported in Ref. 6. 

The unsteady aerodynamic analysis used in this investigation (c.f. Refs. 1 
through 4) applies to fan or compressor cascades operating at subsonic inlet Mach 
numbers. This analysis accounts for the effects of blade geometry, mean blade 
loading and transonic phenomena, including shocks and their motions, on the 
unsteady aerodynamic response to prescribed blade motions. The unsteady equa- 
tions are derived from the assumption that unsteady disturbances are of small- 
amplitude and harmonic in time relative to a fully nonuniform irrotational mean 
or steady background flow. The resulting set of unsteady equations are linear, 
time-independent and contain variable coefficients which depend on the underlying 
mean flow. These equations are solved using an implicit least-squares finite- 
difference approximation which is applicable on arbitrary grids — an important 
feature for turbomachinery applications. In previous work numerical solutions, 
based on this linearized analysis, have been reported for subsonic flows through 
vibrating cascades of double-circular-arc (DCA) airfoils and NACA 0012 airfoils 
(Refs. 1 and 2), and for subsonic and transonic flows through vibrating cascades 
of flat-bottomed DCA airfoils (Refs. 3 and 4). More recently the unsteady 
analysis has been applied to representative two-dimensional outer-span sections 
of an actual fan rotor (Ref. 7). The results reported in Ref. 7 reveal dramatic 
effects of mean blade loading on the unsteady response to thin blades operating 
at high subsonic inlet Mach number and vibrating at high frequency. 

The present study has been motivated by the need to assess the foregoing 
linearized unsteady aerodynamic analysis through comparisons between theoretical 
predictions and available experimental data. This exercise should serve to guide 
the improvements in numerical modeling that will be required to meet the goal of 
providing an efficient and reliable unsteady aerodynamic analysis, which can be 
used in turbomachinery aeroelastic design investigations. 
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THE LINEARIZED UNSTEADY AERODYNAMIC ANALYSIS 


Governing Equations 

We consider isentropic and irrotational flow of a perfect gas through a two- 
dimensional cascade (see Fig. l) of vibrating airfoils. The blades are 
undergoing identical harmonic motions at frequency U), but with a constant phase 
angle o between the motions of adjacent blades. It is assumed that the flow 
remains attached to the blade surfaces and that the blade motion is the only 
source of unsteady excitation. 

As a result of the foregoing assumptions the flow through the cascade is 
governed by the field equations 


3 p ~ ~ 

+ V • (p V $) = 0 

3 t 


(1) 


and 


(p/pj) = 1- (Y-l) Pi {* t + [(vi) 2 - V^]/2}/(TP 1 ) . (2) 


Here ^(X,t) and p(X,t) are the time-dependent velocity potential and fluid 
density; P^ , P]_ and are the upstream free-stream density, pressure and 

velocity respectively, Y is the specific heat ratio of the fluid, it is a position 
vector and t is time. In addition to Eqs. (1) and (2), the flow must be 
tangential to the moving blade surfaces and acoustic energy must either attenuate 
or propagate away from or parallel to the blade row in the far field. Finally, 
we also require that mass and tangential momentum be conserved across shocks and 
that pressure and the normal component of the fluid velocity be continuous across 
the vortex-sheet unsteady wakes which eminate from the blade trailing edges and 
extend downstream. 

Equations (1) and (2) along with the equations based on the foregoing 
conditions at blade, shock and wake surfaces and in the far field are sufficient 
to determine the unsteady flow. However, the computing resources required by 
this nonlinear time-dependent unsteady aerodynamic formulation limits it 
usefulness for turbomachinery aeroelastic investigations. Instead, a small- 
unsteady-disturbance assumption is usually invoked. This assumption permits an 
efficient approximate description of the unsteady flow which is suitable for 
aeroelastic calculations. 
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Thus, following the approach used in. Refs. 1-4, we assume that the blades are 
undergoing small-amplitude (i.e., of £?(e) « 1) unsteady motions and expand the 
flow variables in asymptotic series in £; e.g., 


$ (X,t) = $ (X) + <J> ( x , t ) +..... (3) 

Here, $>(X) is the zeroth-order or steady-flow potential, <j>(^,t) = Re {4) ) e x0Jt } 
is the first-order (in e) unsteady perturbation potential produced by harmonic 
blade motions, the dots refer to the higher-order terms and Re { } denotes the 
real part of { }. In addition to (3), Taylor series expansions are used to refer 
information on moving blade, shock and wake surfaces to the respective mean 
positions of these surfaces. After substituting these expansions into the full 
governing equations, equating terms of like power in £, and neglecting terms of 
higher than first order in e, time-independent nonlinear and linear variable- 
coefficient boundary-value problems are obtained, respectively, for the zeroth- 
and first-order flows. 

The field equations governing the steady flow follow from Eqs. (1) and (2) 
after replacing the time-dependent variables $(S,t) and p(£,t) by their zeroth- 
order or steady-flow counterparts $(X) and p(i) , and setting temporal derivative 
terms equal to zero. The resulting equations, when combined with the associated 
zeroth-order boundary condition of flow tangency at the mean blade surfaces, 
prescribed uniform flow conditions at the inflow boundary, and a Kutta condition 
at blade trailing edges, describe the steady background flow through the 
stationary cascade. 

The differential equation governing the first-order or linearized unsteady 
flow, i.e., 


a 2 V 2 <}> = 



(7-1) V 2 $ 


<f> + V(V4>) 2 * V<t>/2, 
Dt 


(4) 


follows from the mass conservation law (1), Bernoulli’s equation (2), the 
isentropic relations, and the asymptotic expansions for the flow variables. Here 
A - Y P/P is the speed of sound propagation in the steady background flow, 

Dg /Dt = iu) + V$*V is a mean flow convection derivative operator and is the 
complex amplitude of the linearized unsteady potential. Solutions to Eq. (4) are 



subject to both boundary conditions at the mean positions of the blade, shock and 
wake surfaces and requirements on the behavior of the unsteady disturbances far 
upstream and downstream from the blade row. Shock and wake (i.e,, the steady 
downstream stagnation streamlines) mean positions are determined from the steady 
solution. The unsteady surface and far-field conditions are given explicitly in 
Refs. 1~4 and will not be repeated here. 

The preceding aerodynamic formulation for determining the unsteady flow 
through a cascade of airfoils undergoing small-ampl itude harmonic oscillations 
requires the solution of a nonlinear boundary-value problem for the zeroth-order 
or steady flow, followed by the solution of a linear variable-coefficient 
boundary-value problem for the first-order or linearized unsteady flow. Both of 
these problems are time-independent. Moreover, because of the cascade geometry 
and the assumed form of the blade motion, the steady and linearized unsteady 
flows must exhibit blade-to-bl ade periodicity. Thus; for example, 


P(X) . (5) 

and 


P(X + mle n ) = 


p(X + mTe^) = P(X) e im0 , (6) 


where P(J) is the steady pressure, p(^) is the complex amplitude of the 
linearized unsteady pressure, m = 0 , +_ 1 , +_ 2 , ... is a blade number index, T is 

the blade spacing and e^ is a unit vector in the "c ircumf erent ial" or p-direction 
(see Fig. 1) . 

Conditions (5) and (6) allow a numerical resolution of the steady and the 
linearized unsteady flow equations to be restricted to a single extended blade- 
passage region of the cascade. Although the unsteady solution is dependent on 
the steady solution, the numerical procedures used to solve the two equation sets 
can be independent of each other. The numerical approximations used in the 
present study are described briefly below with primary emphasis placed on the 
unsteady solution procedure. 
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Numerical Procedures 


In view of the stringent and often conflicting requirements placed on the 
construction of a computational mesh suitable for the resolution of cascade flow 
fields, numerical approximations to the steady (Refs. 8 and 9) and the 
linearized unsteady (Refs. 1-4) problems have been based on a two-step solution 
procedure. First, large-scale phenomena are determined on a sheared H-type 
cascade mesh of moderate density. Then, for blades with rounded leading edges or 
flows containing shocks, the second step is to determine detailed local solutions 
on body-fitted polar-type meshes of high density. The local mesh domains are 
chosen to cover and extend well beyond limited regions of high mean-velocity 
gradient. The final solutions to the steady and unsteady boundary-value problems 
are taken to be composites of the corresponding cascade- and local-mesh 
solut ions . 

The steady or mean flow solutions reported herein were determined using the 
two-dimensional finite-area approximation described in Refs. 8 and 9. In this 
method the integral form of the continuity equation is approximated over polygons 
in the physical plane. These polygons are constructed by a tr iangularizat ion of 
the computational mesh. For transonic flows artificial compressibility is 
introduced in supersonic regions to stabilize the solution scheme and to capture 
shocks . 

The unsteady solutions were determined using the finite-difference 
approximation described in Refs. 1-4. In this approach, discrete 
approximations to the linear unsteady equations are obtained using an implicit 
least-squares interpolation procedure. Thus, an algebraic approximation, L<|>, to 
the linear differential operator, at the mesh point, Qq, is written in terms 

of the values of <j> at Qq and at certain neighboring points Q^, . . . , Q N as 
follows 


= = Ao + 


N 

l 


n=l 


6 

n 


( +n" V- 


(7) 


The difference coefficients 6^ in Eq. (7) are evaluated in terms of a prescribed 
set of interpolating functions and a set of interpolating coefficients. The 
latter are determined by a weighted least-squares procedure. 

The points Qq through Q^, termed a neighbor set, are defined in a "centered" 
fashion for interior or field points and in a one-sided fashion for boundary 
points. For transonic applications one must distinguish between regions of 
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subsonic flow, where the unsteady differential equation is elliptic, and super- 
sonic flow, where it is hyperbolic. This change in equation type depends on the 
local mean-flow Mach number and is accommodated through the use of local type- 
dependent differencing approximations (see Refs. 3 and 4). With proper ordering 
the discrete unsteady equations can be assembled into a single block- 
pentad iagonal system, which can be solved conveniently using Gaussian 
elimination . 
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UNSTEADY EXCITATION AND RESPONSE PARAMETERS 


In each of the aeroelastic cases considered here the unsteady excitation is 
due to blades undergoing prescribed s ingle-degree-of-f reedom torsional (pitching) 
motions about an axis at or near midchord. Furthermore, the blades are assumed 
to vibrate at a frequency w (or a reduced frequency k based on blade semi-chord) 
with a constant phase angle o between the motions of adjacent blades. Thus the 
angular displacement, ot^ ( t ) , of the mth blade is given by 


* m (t) = Re{a m e lUJt } - 'Re{a Q e 


i(oi)t+ma) 1 


m=0,±l,±2,. 


( 8 ) 


where a m is the complex amplitude of the mth blade angular displacement. 

The linearized unsteady response parameters obey relations similar to (8) 
and thus it is only necessary to consider the response parameters associated with 
the reference blade. For convenience the subscript m = 0, which refers to the 
reference blade, will be omitted in the following discussion. The unsteady 
aerodynamic response quantities of interest include the perturbation unsteady 
pressure and pressure-difference coefficients, C and AC , respectively; the 
unsteady aerodynamic moment coefficient, C , acting an the moving reference blade 
surface; the aerodynamic work done by the airstream on this blade over one cycle 
of its motion, C w ; and the aerodynamic damping, coefficient, S. These 
quantities are defined below. 

The perturbation unsteady pressure coefficient at the moving reference blade 
surface is given by 

C p (x,t) = Re{c p (x)e i “ t } = Re{ |C p j e l(U>t+ V } 

= (x,t) y/ P^ 2 |a|^ , 0 _< x £ c, (9) 

where C^(x) is the complex amplitude of this pressure coefficient, x is a 
coordinate measuring distance along the mean blade-chord line, | | denotes the 

magnitude of a complex quantity and <t> p is the phase angle by which the complex 
vector C (x) leads the complex-displacement vector a. Furthermore, p / ^ ? (x,t) is 
the unsteady pressure perturbation at the instantaneous position ot the reference 
blade surface^, which is determined theoretically from the solutions for the 
steady and linearized unsteady veloicty potentials; and are the uniform 
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density and flow speed, respectively, far-upstream of the cascade (see Fig. 1); 
and c is the blade chord length. The unsteady pressure-difference coefficient is 
defined by 


ACp(x , t) = Re{AC p (x)e ia)t } = Re { j AC p j e l(a)t + <f>Ap) } 


C p( x _ 5 t) “ C (x + ,t) , 0 £ x £ c 


( 10 ) 


where the subscripts - and + refer to the lower (pressure) and upper (suction) 
surfaces of the blade respectively, and 4^ is the phase angle by which the 
complex pressure difference, ACp, leads the complex angular displacement, a. 

The perturbat ion unsteady moment coefficient is defined by 


C (t) 
m 


- 


103 1 ' 




C p (x , t ) (R p *ds ) , 


( 11 ) 


where is the phase angle by which the moment leads the angular displace- 
ment, is a position vector extending from the reference blade axis of rotation 
to a point on the moving reference blade surface and ds is a differential vector 
tangent to this blade surface and directed counterclockwise. Both the angular 
displacement and the moment are regarded here as being positive in the clockwise 
direction. It should be noted that, if shock discontinuities are present, 
additional terms must be added to the right-hand-sides of Eqs. (9) and (11) to 
account for the concentrated loads produced by the motions of the shocks along 
the blade surface (see Refs. 3 and 4). 

Finally, the aerodynamic work per cycle and the aerodynamic damping 
coefficient for the pure torsional vibrations are given by 


C w = — 7T j a j H = Tija 


jcjsin <f> . 
» m 1 m 


( 12 ) 


When a is prescribed as a real quantity, the right -hand-s ide of Eq. (12) can be 
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replaced by 7ralm{c m }, where Im{ } denotes the imaginary part of { } . The 
stability of the torsional blade motion (according to linearized theory) depends 
upon whether C w y- 0 (or H 0) . Thus, if C w < 0, the airstream removes energy 
from the blade motion and this motion is stable; if = 0, there is no net 
transfer of energy and the blade motion is neutrally stable; and finally, if 
> 0, the airstream supplies energy to the blade motion and this motion is 
unstable . 
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THE STANDARD CASCADE CONFIGURATIONS 


At the Second International Symposium on Aeroelasticity in Turbomachines 
held in Lausanne, Switzerland in September 1980 (Ref. 10), it was proposed that 
an experimental and analytical data base be established for selected turbo- 
machine-cascade configurations. The purposes of this effort are: first, to 

provide a convenient means for evaluating and comparing the results of different 
theoretical analyses; and second, to provide a convenient means for assessing 
theoretical predictions through comparisons with experimental measurements. In 
the four year period between the Second Symposium and the Third Symposium, held 
in Cambridge, U.K. , nine two-dimensional and quasi three-dimensional standard 
cascade configurations were proposed, and the developers of various theoretical 
models were asked to apply their calculation methods to predict the aeroelastic 
behavior of these configurations. 

The present authors have contributed to this project by applying the un- 
steady aerodynamic analysis of Refs. 1 through 4 to predict unsteady aerodynamic 
response properties for the following two-dimensional cascade configurations: 

(1) a compressor cascade consisting of cambered NACA 65 series airfoils and 
referred to in Ref. 5 as the first standard configuration; 

(2) a subsonic/transonic fan-tip cascade, referred to as the fifth standard 
conf igur at ion; 

(3) a cascade of flat plate airfoils, referred to as the eighth standard 
configuration; and 

(4) a cascade of double-circular-arc airfoils referred to as the ninth 
standard configuration. 

These standard configurat ions^ along with the recommended aeroelastic test cases, 
are described in detail in Ref. 5. We determined unsteady aerodynamic response 
information for fifteen (15) recommended test cases on the first, eighth and 
ninth configurations, and for eight (8) cases on the fifth configuration for a 
total of fifty-three (53) test cases. 

Our results were submitted to T. H. Fransson for inclusion in a paper (Ref. 
11) that he presented at the Third Symposium. In addition, one of the present 
authors (J. M. Verdon) attended this Symposium to give a brief presentation 
describing the linearized unsteady potential-flow analysis used in our study and 
the application of this analysis to the four standard configurations listed 
above. A sample of our results are included in Figs. 1 through 4 of Ref. 11. 
These figures contain our predictions (referred to as Code 3) as well as the 
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analytical predictions and the experimental data contributed by other 
investigators. A more complete report on the analytical and experimental results 
submitted on the standard configurations is being prepared by Fransson for 
publication later this year. In view of this, a detailed presentation of our 
results will not be given in this report. But we will include a representative 
selection of these results, including several based on more recent calculations, 
and offer some general observations on the analytical and experimental results 
presented at the Third Symposium for the standard configurations listed above. 


The First Standard Configuration 

The first standard configuration is a compressor-type cascade operating at 
low subsonic inlet Mach number (M ^ ) . This cascade has a 35° stagger angle ( 0 in 
Fig. 1) and a gap to chord ratio t/c of 0.73. The blades are constructed by 
superimposing a NACA 65 series thickness distribution on a 10 degree circular-arc 
camber line. They are driven in the torsional mode at amplitudes, |a|, of 0.5° 
and 2.0° around an axis at (x/c,y/c) = (0.5, 0.115) and at frequencies varying 
from 9 to 77 Hz (i.e., semi-chord reduced frequencies, k, between 0.072 and 
0.603). Most of the tests were conducted at 15.5 Hz (k = 0.122). The experi- 
ments were conducted by F. 0. Carta (Ref. 12) in the UTRC Oscillating Cascade 
Wind Tunnel and three investigators contributed theoretical predictions. 

In general, the theoretical analyses predicted the measured unsteady 
aerodynamic response behavior quite well, both in terms of the detailed unsteady 
pressure information along a blade surface and the global unsteady response 
quantities (unsteady moment, etc.). However, there is some ambiguity concerning 
the inlet flow angle (6^). Theoretical predictions for mean surface pressures 
have been found to match closely the corresponding measured pressures, but only 
if the inlet flow angle prescribed for the theoretical calculation differs from 
that reported in the experiment. The present authors have found that a similar 
situation also exists for the fifth standard configuration. 


A sample of our predictions along with Carta* s measurements for the mean 
pressure coefficient, Cp , on the reference blade surface, the unsteady pressure 
coefficient, the unsteady pressure-difference coefficient, and the unsteady 
moment are shown in Figs. 2 through 5, respectively. The mean pressure 
coefficient is defined by 


C p (x) = (p(x) 


- p i)/(i Vi 2 ) 


(13) 


The analytical and experimental results shown in Figs. 2 through 3 are in very 
good agreement. They were determined for an inlet Mach number of 0.17 and blades 
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vibrating with a = 2° and at a reduced frequency of 0.122. The inlet flow angle 
used in the analysis is 27.8°, while that reported for the experiment is 24°. 

It should be noted that the theoretical moment curve shown in Fig. 5 was deter- 
mined recently on the basis of moment predictions for forty values of 0, so that 
a detailed description of the unsteady moment behavior, especially in the super- 
resonant region near 0 = 0°^ could be provided. Results for only nine values of 
0, corresponding to those at which measurements are available, were submitted at 
the Third Symposium. 


The Fifth Standard Configuration 

The fifth standard configuration is a two-dimensional subsonic/ transonic 
fan tip cascade which has been tested in a rectilinear cascade wind tunnel at 
ONERA by Szychenyi and Finas (Ref. 13). The experimental configuration consists 
of six thin uncambered blades having a maximum thickness to chord ratio of 0.027 
at approximately 67% of blade chord. The cascade stagger angle, 9, is 30.7° and 
the gap to chord ratio, t/c, is 0.95. The experiments were performed for 
torsional motions about midchord with |a| = 0.3° and over a frequency range 
extending from 75 (k= 0.14) to 550 (k = 1.02) Hz. The inlet Mach number was 
varied from 0.5 to 1.0, and the inlet flow angle was varied from 28.7° to 18.7°, 
corresponding to attached through fully separated flow. Unfortunately, only one 
blade was vibrated in the experiments. Although methods are available for con- 
verting the resulting unsteady response information to represent that for a cas- 
cade with all blades vibrating, these were not applied and therefore some ques- 
tions remain concerning the validity of the comparisons between the data and the 
analytical predictions for the fifth standard configuration. 

Two theoretical models (c.f. Refs. 1-4 and Ref. 14) were applied to provide 
predictions for = 0.5 in the attached (experimental 6^ = 26.7°) and partly- 
separated (experimental 6^ = 24.7°) flow domains. No predictions were submitted 
for the fully-stalled experimental conditions. There are large differences 
between the predicted and measured values for the magnitudes of the unsteady 
pressures, pressure differences and aerodynamic moments. Moreover, substantial 
differences also exist between the results of the two theoretical prediction 
methods for these quantities. However, the two analytical and the experimental 
results for the phase angles between the unsteady surface pressure, pressure 
difference and aerodynamic moment and the angular displacement are in reasonably 
good agreement. Consequently, both theory and experiment give the same indica- 
tion in regard to the stability of the torsional blade motions. The differences 
between the theoretical predictions and the measurements can perhaps be 
attributed to viscous separation phenomena which are present in the experiment 
but are not accounted for in the theoretical models. In addition, only one blade 
was excited in the finite-cascade experiment while all blades are assumed to be 
oscillating in the infinite-cascade theoretical formulations. The reasons for 
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the differences between the predictions of the two theoretical models are not 
apparent at the present time. 

A selection of our results for mean surface-pressure coefficient, unsteady 
pressure-difference coefficient and aerodynamic damping are shown in Figs. 6 
through 8 along with the data contributed by Szychenyi. Substantial differences 
exist between the predictions and measurements for the mean-pressure coefficient 
near the leading edge of the blade (Fig. 1) and for the magnitude of the unsteady 
pressure-difference coefficient over the entire blade. Good agreement exists 
between the theoretical and experimental values for the phase angle, , and the 
aerodynamic damping, S 3 for the partly-separated flow, but this must be regarded 
as fortuitous in view of the discrepancies between the analytical and 
experimental results for |ACp(x)|. 


The Eighth and Ninth Standard Configurations 

The eighth and ninth configurations were selected for the validation of 
various theoretical prediction methods, especially at high subsonic inlet Mach 
number. Therefore, no experimental results were submitted for these standard 
configurations. The eighth configuration is a two-dimensional cascade of flat- 
plate blades. The blades undergo torsional oscillations about midchord at 
relatively high reduced frequency (i.e., k = 1.0), and the phase angle between the 
motion of adjacent blades is 90°. The cascade stagger angle and gap/chord ratio 
are variable, but for most of the test cases 6 = 30° and t/c = 0.75. The inlet 
Mach number varies from 0 to 1.5, and in each case the inlet flow angle is set 
equal to the cascade stagger angle and therefore the mean flow is uniform. 

The present authors submitted results only for subsonic inlet Mach numbers 
with M| extending from 0.1 to 0.95. Three other investigators also submitted 
theoretical predictions for this Mach number range. The results of the 
various prediction methods for the subsonic flat-plate cascades are generally in 
good agreement. However, there are some discrepancies between the predictions 
for |Cp| at high subsonic inlet Mach number. These are due, perhaps, to the 
difficulties associated with the numerical prediction of high-frequency unsteady 
flow at high subsonic Mach number. Our results for the unsteady pressure- 
difference coefficient at M| - 0.8 and unsteady moment versus inlet Mach number 
are given in Figs. 9 and 10, respectively. The results of the other prediction 
methods are in good agreement with those presented in these figures. 

The ninth configuration is a two-dimensional cascade of double-circular-arc 
airfoils. Both symmetric and flat-bottomed DCA airfoils have been recommended. 

In the symmetric case the blade thickness varies from one to ten percent of 
chord, while the flat-bottomed blades are 5% thick. The cascade gap to chord 
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ratio is 0.75 and there are two stagger angles, 9 = 30° and 9 = 45°. The blades 
undergo torsional oscillations about midchord at a reduced f requency^ of 1.0. 

The inlet Mach number varies from 0 to 1.5, and the inlet flow angle is deter- 
mined by prescribing a unique mean-incidence condition at blade leading edges. 

The results of two theoretical analyses (Refs. 1-4 and Ref. 14) were sub- 
mitted for subsonic flow through the symmetric DCA configurations. The present 
authors also contributed results for subsonic and discontinuous transonic flows 
through the cascade of flat-bottomed DCA airfoils; however, we did not submit 
results for supersonic inlet conditions. The results of the two analyses for 
subsonic flows through the cascades of symmetric DCA blades are generally in good 
agreement, but increasing differences are observed with increasing inlet Mach 
number and/or blade thickness. Our results for the unsteady pressure-difference 
coefficient and for the unsteady moment versus inlet Mach number for a cascade of 
2% thick blades are shown in Figs. 11 and 12. The results of the other predic- 
tion method are in close agreement with those shown in these figures. Again, we 
note that the moment curves shown in Figs. 10 and 12 are based on recent more 
detailed calculations than those submitted for the Third Symposium. 
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THE NASA LEWIS FLUTTER CASCADE 


The NASA Lewis flutter cascade consists of nine uncambered or symmetric 
biconvex airfoils which are driven simultaneously to provide a 1.2° amplitude 
pitching motion about midchord. The cascade has a stagger angle of 37° and a 
gap/chord ratio of 0.767. The symmetric biconvex airfoils are 7.62 cm. long, and 
they have a maximum thickness at midchord of 0.58 cm. The radius of curvature of 
the upper (suction) and lower (pressure) surfaces is 27.4 cm. The experimental 
airfoils close in rounded leading and trailing edge sections of radius 0.000254 
cm. Tests were performed at an inlet Mach number of 0.65 and an inlet flow angle 
of 0° in an attempt to simulate an unloaded and shock-free mean-flow condition. 
Tests were also performed at an inlet Mach number of 0.8 and an inlet flow angle 
of 30° in order to observe the surface-pressure response when an oscillating 
shock occurs near the leading edge of each airfoil. The unsteady response data 
are given in Ref. 6 for torsional vibrations about midchord at interblade phase 
angles of -90° and +90° and vibration frequencies of 200 Hz and 500 Hz (i.e., 
semi-chord reduced frequencies of k = 0.221 and 0.553 for Mj, = 0.65, and k = 

0.183 and 0.458 for M ] _ - 0.80). 

Analytical results for the NASA Lewis Cascade are depicted, along with the 
measurements of Ref. 6, in Figs. 13 through 26 of the present report. The blades 
used in the analytical study are extended versions of those used in the experi- 
ment. That is, sharp-edged biconvex (or symmetric DCA) blades having the same 
surface radius and the same thickness at midchord as those used in the experi- 
ment were employed, but the upper and lower surfaces of the blades used in the 
analysis were extended to close in wedge-shaped leading and trailing edges. As a 
result the blades used in the analysis are 4% longer than those used in the 
experiment. This modification was introduced to simplify analytical considera- 
tions at blade edges. With the exception of the inlet flow angle for the M ^ = 
0.65 case, all other experimental parameters have been retained in the analytical 
study. The steady full-potential analysis (Refs. 8 and 9) used in this study 
indicates that, for = 0.65, the unloaded mean-flow condition occurs at 3 ^ = 
33.22° and not at 3| = 37°, as reported in the experiment. After a number of 
unsuccessful attempts to match the measured blade-surface Mach number distribu- 
tions by prescribing different combinations of inlet and exit flow angles, we 
decided to provide analytical predictions at M^ = 0.65 for the unloaded mean-flow 
condition. These results should provide useful "benchmark" information for use 
in possible future experiments. 

The predicted and measured unsteady response information for M^ = 0.65 are 
shown in Figs. 13 through 19. As mentioned, the analytical predictions corres- 
pond to an inlet flow angle of 33.22° while the inlet flow angle reported in the 


17 



experiment is 37.0°. The predicted steady Mach number distributions in Fig. 13 
correspond to the zero mean blade-loading condition, while the experimental Mach 
number distributions indicate that a negative steady lift force acts on each 
blade. Although the analytical steady background flow differs from that 
occurring in the experiment, the mean-flow gradients are relatively small in both 
instances and are not expected to have a major impact on the unsteady response. 

It is therefore still meaningful to compare the predicted and measured unsteady 
response information. 

The predicted and measured unsteady surface pressures are shown in Figs. 14 
through 17 for the following combinations of frequency and interblade phase 
angle: o> = 200 Hz and o = -90° (Fig. 14); oj = 200 Hz and a = 90° (Fig. 13); w = 
500 Hz and a = -90° (Fig. 16); and o> = 500 Hz and O = 90° (Fig. 17). With the 
exception of the latter combination, which is close to a resonance condition (see 
Fig. 19), a qualitative agreement between the measured and predicted unsteady 
surface pressures has been achieved. The predicted behavior of the unsteady 
moment coefficient versus interblade phase angle is illustrated, along with the 
experimental values at 0 = ±90°, in Fig. 18 for 0) = 200 Hz and in Fig. 19 for 
0 ) = 500 Hz. It should be noted that the experimental moment coefficients were 
estimated from unsteady pressure-response data available at only six points on 
the blade suction surface and six points on the blade pressure surface. Hence, 
the experimental estimates for the unsteady moment coefficients may not be too 
reliable. However, with the exception of the motion at to = 200 Hz and 
0 = 9o° (Fig. 18), the in-phase components (Re{C m }) of the analytical and experi- 
mental moments are in good agreement. The out-of-phase components ■( Im{C m }) are 
also in good agreement for o = 90°, but appreciable differences exist at a ~ -90° 
(see Figs. 18 and 19). 

Similar results are provided in Figs. 20 through 26 for the high inlet Mach 
number (M^ - 0.8), high mean incidence (6^ = 30°) condition. In this case the 
analytical predictions in the vicinity of the blade leading edge are questionable 
because the steady and unsteady flow behavior in the vicinity of a sharp leading 
edge at incidence cannot be predicted adequately. Indeed, unless a relatively 
coarse computational grid is employed near the leading edge of a blade, iterative 
solution procedures for the steady full-potential equation will diverge. It is 
therefore not possible to predict the unsteady effects associated with certain 
phenomena observed in the experiment (Ref. 6), e.g., leading-edge lambda-type 
shocks and flow separations downstream of these shocks, using the current 
inviscid steady and unsteady codes. 

In view of these limitations the correlation between the predicted and 
measured surface Mach number distributions, shown in Fig. 20, is encouraging. 

The predicted and measured unsteady surface-pressure coefficients for torsional 
vibrations about midchord, at the same combinations of frequency and interblade 
phase angle considered above, are shown in Figs. 21 through 24. Again, with the 
exception of the motion at to = 500 Hz and a = 90°, which is near a resonance 
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condition, there is a reasonable qualitative agreement between the predicted and 
measured surface-pressure coef f ic lent s . However, large differences do exist 
over the forward part of the blade suction surface. The predicted and measured 
(at 0 = -90°) unsteady moment coefficients for the flows at - 0.8 are shown 
in Fig. 25 for u> = 200 Hz and in Fig. 26 for a) = 500 Hz. 
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CONCLUDING REMARKS 


This report describes a theoretical contribution to the establishment of a 
theoretical and experimental data base for the validation of current and future 
unsteady flow analyses. In particular, the linearized potential-flow analysis of 
Refs. 1 through 4 has been applied to provide unsteady aerodynamic response 
predictions for five two-dimensional compressor (or fan) - type cascade configu- 
rations. These include the first, fifth, eighth and ninth standard cascade 
configurations of Ref. 5 and the NASA Lewis flutter cascade of Ref. 6. Unsteady 
response measurements are available for the first and fifth standard configura- 
tions and for the NASA Lewis cascade; several other investigators have contributed 
theoretical predictions for the standard cascade configurations. 

Our predictions have been found to be in very good agreement with those of 
other theoretical investigators for the first (a cascade of NACA 65 series 
airfoils), eighth (a flat-plate cascade) and ninth (a cascade of double-circular- 
arc airfoils) standard configurations (see Ref. 11), and with the experimental 
measurements of Ref. 12 for the first standard configuration. In addition, good 
qualitative agreement exists between our predictions and the experimental 
measurements of Ref. 6 for the NASA Lewis configuration (a cascade of symmetric 
biconvex airfoils). However, wide discrepancies exist between the predicted and 
measured values for the magnitudes of the complex response quantities for the 
fifth standard configuration. These are due, possibly, to the different 
conditions which prevail in the experimental (only one blade vibrating) and the 
analytical (all blades vibrating) flow environments. 

In both the NASA Lewis and the fifth standard configurations the blades have 
very small leading-edge radii, and they operate at non-zero mean incidence. 

Thus, large flow gradients occur along each blade suction surface in the vicinity 
of its leading edge, and these promote shock formation and flow separation. At 
present it is not possible to resolve such detailed local phenomena with existing 
inviscid (or viscous) flow codes. Thus, we recommend that future research be 
directed toward improving inviscid steady and linearized unsteady flow codes so 
that transonic leading-edge phenomena can be predicted. In addition, the 
possibility of extending linearized unsteady aerodynamic formulations to account 
for viscid/ inviscid interaction phenomena, including leading-edge and/or shock- 
induced separations, should be explored. 

The concept of establishing a theoretical and experimental data base for 
standard cascade configurations, as proposed at the Second International 
Symposium on Aeroe las t ic i ty in Turbomachines, is an important one. Before the 
start of the Third Symposium, there were concentrated efforts on the part of both 
experimentalists and theoreticians, working independently, to submit measurements 
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and predictions for the various configurations. Such an approach is useful in 
getting started, especially if the measurements and predictions are found to be 
in close agreement. Usually this will not be the case, however, because of the 
complicated nature of the unsteady flows through turbomachinery blade rows. 

Thus, future activities of this kind would benefit from research programs in 
which theoreticians and experimentalists worked together from the outset. In 
such programs, benchmark analytical and experimental studies could be planned and 
conducted to identify carefully both the regions in which the results of theory 
and experiment will correlate and the additional physical phenomena that must be 
incorporated in theoretical models. Basic issues to be addressed are the deter- 
mination of the ranges of validity of both the linearized unsteady aerodynamic 
analyses which have been developed for turbomachinery aeroelastic applications 
and the finite-cascade wind-tunnel experiments which are used to validate these 
analyses . 
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Speed of sound propagation in the steady background flow, Eq. (4) 
Perturbation unsteady moment coefficient, Eq. (11) 

Perturbation unsteady surface pressure coefficient, Eq. (9) 

Steady pressure coefficient, Eq. (13) 

Complex amplitude of the perturbation unsteady moment coefficient, 
Eq. (11) 

Complex amplitude of the perturbation unsteady surface pressure 
coefficient, Eq. (9) 

Aerodynamic work per cycle of blade motion, Eq. (12) 

Blade chord length. Fig. 1 

Mean-flow convective derivative operator, Eq. (4) 

Differential vector tangent to blade surface and directed 
counterclockwise, Eq. (11) 

Unit vector in the cascade "c ircumf erent ial" or redirection, Eq. (5) 
Imaginary part of { } 

Reduced frequency based on blade semi-chord 
Linear differential operator, Eq. (7) 

Linear difference operator, Eq. (7) 

Zeroth-order or steady Mach number 

Blade number index (m=0 denotes the reference blade) 

Number of neighbor mesh points, Eq. (7) 

Mesh point index, Eq. (7) 

Zeroth-order or steady pressure, Eq. (3) 
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Linearized unsteady pressure at the instantaneous position of the 
reference blade surface, Eq . (9) 

Complex amplitude of first order or 1 inearized uns teady pressure, 
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Mesh point, n=0 refers to calculation point, n=l,..., N refers to 
neighboring mesh points 

Multiplicative constant , Eq . (7 ) 

Position vector extending from reference blade axis of rotation to 
points on the reference blade surface, Eq . (11) 

Real part of { }, Eq. (3) 

Time 

Zeroth-order or steady velocity , Eq. (2) 

Position vector 

Distance from leading edge along blade chord, Fig. 1 
Blade angular displacement, Eq. (8) 

Complex amplitude of the blade angular displacement, Eq. (8) 

Inlet, exit flow angle, Fig. 1 
Specific heat ratio, Eq. (2) 

Perturbation unsteady pressure-difference coefficient, Eq. (10) 

Complex amplitude of the perturbation unsteady pressure-difference 
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Finite-difference coefficient, Eq. (7) 

Small parameter which is a measure of the amplitude of the blade 
motion 

Cascade "circumferential” coordinate, Fig. 1 
Cascade stagger angle, Fig. 1 
Aerodynamic damping coefficient, Eq. (12) 

Nonlinear time-dependent fluid density, Eq. (1) 

Zeroth-order or steady density, Eq. (1) 
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Nonlinear time-dependent velocity potential, Eq . (1) 
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to Frequency of the blade motion, Eq. (8) 

Subscripts 


_B Instantaneous position of reference blade 

m Blade number index, Eq. (8) 
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n Mesh point index, Eq. (7) 

t Partial derivative with respect to time, Eq.(l) 

1 Uniform inlet flow condition, Eq. (2) 
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ANALYSIS: 3 1 =27.8° 

UPPER SURFACE 

LOWER SURFACE 


EXPERIMENT: /3 1 =24° 
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Fig. 2 Mean Surface-Pressure Coefficient for the First Standard 
Configuration; M 1 =0.17. 



















ANALYSIS: ^ =28.2° EXPERIMENT: ^ =26.7° 
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Fig. 6 Mean Surface-Pressure Coefficient for the Fifth Standard 
Configuration; M-j =0.5. 















Fig. 10 Unsteady Moment Coefficient Versus inlet Mach Number for 
the Eighth Standard Configuration; /3-j = 0 = 30°, k = 1.0, ct = 90 












Fig. 12 Unsteady Moment Coefficient Versus Inlet Mach Number 
for the Ninth Standard Configuration (2% Thick 
Symmetric Blades); k = 1.0, a = 90°. 












ANALYSIS'. j3 1 =33.22° EXPERIMENT: 0 ^= 37 ° 
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Fig. 15 Unsteady Surface-Pressure Coefficient for the NASA 
Lewis Cascade; M-| =0.65, w = 200 Hz, <7 = 90°. 






















ANALYSIS 
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Fig. 20 Steady Surface Mach Number Distributions for the 
NASA Lewis Cascade; M 1 =0.8, /S-j =30°. 






















Fig. 25 Unsteady Moment Due to Torsion About Midchord for the 
NASA Lewis Cascade; M-| =0.8, /3-j =30°, a; = 200 Hz. 
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